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Abstract 

We have searched for neutrinoless r lepton decays into £hh or £V^, where £ stands 
for an electron or muon, h for a charged hght hadron, tt or K, and for a neutral 
vector meson, p", i^*(892)'' and cj), using a 158 fb~^ data sample collected with the 
Belle detector at the KEKB e+e~ collider. Since the number of events observed 
are consistent with the expected background, we set upper limits on the branching 
fractions in the range of (1.6 — 8.0) x 10~^ for various decay modes at the 90% 
confidence level. 

Key words: TAU Lepton Flavor Violation 
PACS: 11.30.-j, 12.60.-i, 13.35.Dx, 14.60.Fg 



1 Introduction 

In the Standard Model (SM), lepton-flavor- violating (LFV) decays of charged 
leptons are forbidden, or highly suppressed even if the effect of neutrino mix- 
ing is taken into account [1]. In contrast, LFV decay processes are expected 
to appear with much larger branching fractions than those in the SM if there 
are contributions from new physics. Searches for LFV decay processes may 
thus reveal new physics beyond the SM. Some models predict LFV decays 
of r leptons at a level accessible at the high luminosity S-factories [2,3]. In 
this paper, we report on a search for LFV in fourteen r~ decay modes into 
neutrinoless final states with one charged lepton i and two charged pseu- 
doscalar mesons h: e~n~^7i~ , e~^n~7i~ , n~n~^n~ , /i+7r~7r~, e~Ti'^K~' , e'^n~K~^, 
e+n-R-, e~K+K-, e+R-R-, i^-n+R-, ix-tt-R+, /i+n-R-, /I'R+R- and 
lji'^R~R~ , and eight modes in which r decays into one lepton and one vector 
meson: g-p^, e-R*{892)°, e-R*{892f, (Td, fjrp^, /x-/C*(892)°, i2-R*{892)^ 
and fi^(f). ^ Current upper bounds on the branching fractions for these decays 
are of the order of 10^® at 90% confidence level (CL) and have been set in the 
CLEO experiment using a data sample of 4.79 fb~^ [4]. Very recently CLEO 
results on T — > £hh modes were improved by the BaBar experiment and upper 

^ Charge conjugate decay modes are implied throughout the paper. 
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limits in the range (0.7 — 4.8) x 10~^ were obtained from a 221.4 fb~^ data 
sample [5]. We present here results of a new search based on a data sample 
of 158.0 fb~^ corresponding to 140.9 x 10^ r-pairs collected with the Belle 
detector [6] at the KEKB asymmetric-energy e^e~ coUider [7] operating at or 
near the T(45') resonance. 



2 Event Selection 

The Belle detector is a general purpose detector with excellent capabilities for 
precise vertex determination and particle identification. Tracking of charged 
particles is performed using a three-layer double-sided silicon vertex detec- 
tor (SVD) and a fifty-layer cylindrical drift chamber (CDC) located in a 1.5 
T magnetic field. Charged hadrons are identified by combining dE/dx in- 
formation from the CDC, signal pulse-heights from aerogel Cerenkov coun- 
ters (ACC) and timing information from timc-of-fiight scintillation counters 
(TOF). Photons are reconstructed using a Csl(Tl) electromagnetic calorime- 
ter (ECL). Muons are detected by fourteen layers of resistive plate counters 
interleaved with iron plates (KLM). 

We use TAUOLA [8] for Monte Carlo (MC) event generation of r-pair signals 
and KKMC [9] to implement initial and final state radiation. The KKMC MC 
program predicts a cross-section oia{e^e~ — > r^r~) = (0.8916±0.0006) nb at 
the center-of-mass energy of KEKB. We calculate the number of r-pair events 
from the cross section and measured integrated luminosity £int = 158.01 fb~^. 
The MC data is processed through the Belle detector simulation program 
based on GEANT3 [10] to determine signal efficiencies. We use the CLEO 
QQ event generator [11] for hadronic events and AAFHB [12] for two-photon 
events, and study their contributions to the background of each r decay mode. 

We search for r-pair events in which one r decays into the Ihh (3-prong) final 
state. The other r dominantly decays into one charged particle and any number 
of neutrals (1-prong) with a branching fraction of 85.35% [13]. We require that 
there be four charged tracks with zero net charge and any number of photons 
in an event. We reconstruct the trajectory of a charged track from hits in the 
SVD and CDC, and require that a reconstructed transverse momentum be 
larger than 0.1 GcV/c and polar angle 9 be within the range 25° < 9 < 140°, 
with respect to the direction opposite to the beam. For all charged tracks, 
the distance of the closest approach to the interaction point (IP) is required 
to be within 1 cm transversely and 3 cm along the e+ beam. Photons are 
selected from neutral ECL clusters with an energy threshold E'ciustcr > 0.1 
GeV and are separated by at least 30 cm from the extrapolated projection 
point of any charged track. The tracks and photons in an event are divided 
into two hemispheres in the e'^e~ center-of-mass system (CMS), with a plane 
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perpendicular to the thrust axis calculated from the momenta of all charged 
tracks and photons. We select 3-prong vs. 1-prong topology events, i.e. three 
charged tracks are in one hemisphere and one charged track in the other. We 
define the former hemisphere as the signal side and the latter as the tag side. 
The number of photons on the signal side should be less than or equal to two, 
to allow for photons from initial and final state radiation or photons radiated 
from electron tracks. 

Electrons are identified by an electron likehhood that includes the dE/dx 
value measured with the CDC, the ratio of the chister energy from the ECL 
to the track momentum measured with CDC, ACC hits and shower shape 
in ECL [14]. With this likelihood, electrons are selected with an average effi- 
ciency of 85% over the whole momentum range of those in the signal MC. The 
momentum of a charged track in the laboratory system p should be greater 
than 0.3 GeV/c to avoid poorly identified electrons. In order to correct for 
the energy loss due to bremsstrahlung, the momentum of an electron track 
is recalculated adding the momentum of radiated photon clusters if an ECL 
cluster with energy less than LO GeV is detected within a cone angle of 10° 
along the electron fiight direction. 

The muon likelihood is formed from two variables: the difference between 
the range calculated from the momentum of the particle and that measured 
with the KLM, as well as the value of the KLM hits with respect to the 
extrapolated track [15]. The average efficiency is evaluated as 90% from signal 
MC for muon in the momentum range for LEV decays. For muons, p should 
be larger than 0.6 GeV/c for the same reason as in the electron case. 

Tracks that do not satisfy the requirements for electron or muon candidates 
are classified as hadrons. To distinguish kaons from pions, we use a likelihood 
ratio which is calculated from dE/dx, time-of-flight and the hits in the ACC. 
We achieve kaon efficiencies of 85% in the barrel and 80% in the endcap region. 
In these cases, p should be larger than 0.5 GeV/c. The remaining tracks are 
treated as pions. 

Vector mesons are reconstructed in the following decay modes: p° —>■ tt+tt", 
i^*(892)° K^TT~ and K^K~ . We calculate the vector meson mass. 
My. We then fit the My distribution of the signal MC with two Gaussian 
distributions to take into account the effects of the intrinsic width of the 
resonances and the detector resolution. 

The width of the signal windows for each decay modes is ±1.64(7, where a is the 
standard deviation of the broader Gaussian component: 445 MeV/c^ < Mp < 
1092 MeV/c2, 730 MeV/c^ < Mk* < 1064 MeV/c^ and 1005 MeV/c^ < < 
1035 MeV/c2. 

The background that remains after applying all the selection criteria based on 
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event topology and particle identification is dominated by Bhabha events. This 
background is effectively reduced by requiring the invariant masses calculated 
for all combinations of two oppositely charged tracks to be greater than 0.2 
GeV/c^, assuming the electron mass. 

Backgrounds from BB and cc events are suppressed by the selection criteria 
listed above, however the background from e~^e~ qq continuum events with 
the light quarks u, d and s {uds continuum) is not negligible. Since these 
backgrounds often include 7r° decays, we require that the number of photons 
on the tag side n-y not exceed 1. If the tag is an electron or muon, the 
requirement effectively suppresses the uds continuum background. For tags 
where the tag side is a pion, the selected events are investigated further as 
follows. 

We measure the r flight length l^. using information on IP and r vertex position 
that is reconstructed from the tracks on the signal side. Since r leptons travel 
significant distance before decaying (cr = 87/im), this feature can be used 
to suppress background. Since r-pair events have a more jet-like shape than 
uds continuum events and are distinguishable by use of shape variables. We 
use a normalized second Fox- Wolfram moment R2 [16] to represent the event 
shape. Figure 1 shows the and R2 distributions for r-pair and uds continuum 
events. We calculate the two-dimensional probability density function (PDF) 
in the l^- vs. R2 plane, and form a likelihood ratio Crr/iJ^TT + ^uds), where 
Ct-t and Cuds are the PDFs for r-pair signal and uds continuum background, 
respectively. The condition Ct-t/ {Ct-t + ^uds) > 0.45 is optimized using MC. 
This condition removes 60% of uds continuum background, while 90% of the 
signal is kept when the track on the tag side is a pion. After combination 
of all selections including the and likelihood ratio requirements, the uds 
continuum background is suppressed by a factor of 10^. 

We calculate the missing four-momentum for each event. In signal events, it is 
possible to reconstruct a r mass from the missing momentum, 1-prong charged 
track momentum and all momenta of photons on the tag side, since there is 
no neutrino emission from the signal decays. Figure 2 shows the distributions 
of the r mass on the tag side Mipr for the signal events and uds continuum 
background MC. There is a clear peak at the r mass in the signal distribution, 
while the distribution for background process is smooth. Comparing the yield 
of these distributions, we optimize the central value and width of the signal 
windows independently for electron and muon modes: 0.85 GeV/c^ < Mipr < 
2.31 GeV/c^ for the r ehh and r eV° modes, and 1.06 GeV/c^ < Mipr < 
2.22 GeV/c^ for the r nhh and r nV^ modes. 

To identify signal r decays, we reconstruct the invariant mass M^hh and CMS 
energy Ely^y^ of selected events. We then calculate the energy difference A£^* = 
Elf^f^ — -Ebeam and mass difference AM = M^^^ — M,-, where -Ebeam is the 
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Fig. 1. (a) Flight length It and (b) R2 distributions of the signal MC events (thick 
line) and uds continuum events (thin line) . The histogram normalization is arbitrary. 

beam CMS energy and is the nominal r mass. The signal events should 
have AE* and AM around zero. We define the signal region in the AE*-AM 
plane using the signal MC. Signal MC distributions have tails on the lower 
sides of AE* and AM due to effects from initial and final state radiation, 
bremsstrahlung of electrons and ECL energy leakage. We consider signal MC 
events as properly reconstructed if they satisfy the conditions —0.68 GeV 
< AE* < 0.32 GeV and -0.25 GeV/c^ < AM < 0.25 GeV/c^. We thus define 
a rectangular signal region that contains 90% of properly reconstructed signal 
events. The boundaries of the signal region are summarized in Table 1. 



3 Results 

The signal MC events are generated assuming a uniform decay angular distri- 
bution of r. Signal detection efficiencies e, evaluated from the MC, are listed 
in the fourth column of Table 2 and vary from 2.68% to 5.30%. The actual 
decay angle distribution, however, depends on the model for the LEV interac- 
tion. In order to evaluate the possible effect of correlations, we examine V — A 
and V + A interactions using the formulae given in Ref. [17] and the relative 
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Fig. 2. Mipr distributions for the signal (thick Hne) and uds continuum (thin hne) 

MC events. The left figure (a) is for modes in which an electron is on the signal side, 
and the figure on the right (b) is for muon tags. Dashed lines show the boundaries 
of the signal region. 

differences in the efficiencies from the uniform distribution (Ae/e) are taken 
as systematic errors of detection efficiencies. 



After all selection requirements, some events remain in the signal region. Prom 
MC, we evaluate the number of background events from uds continuum and 
r-pair production. We then scale to the data by the factor aCi^^^/NyiQ, where 
a is the cross section of this process, and A^mc is the number of generated 
events of the process. In the t~ — > e~h'^h'~ and t" — > e'V^ modes, there 
is a contribution from two-photon e+e" — > e'*'e~e'^e~, e^e~ e+e"/^"'"/^", 
e+e" — s> e^e~uu, e~^e~ e^e^dd and e^e~ e^e~ss processes. Since the 
equivalent luminosity for the MC simulation of two-photon processes is much 
smaller than i2int, the contribution of the two-photon processes is estimated 
from a fit to the data in the AM sideband region, in which the contribu- 
tions of uds continuum and r-pair are fixed from MC and the two-photon 
shape is taken from MC with the normalization floated. For all modes the 
sideband region is —0.5 GcV/c^ < AM < 0.5 GcV/c^, while the signal region 
-96 MeV/c^ < AM < 30 MeV/c^ is blinded. The shape of the AM distribu- 
tion is obtained after applying only the particle identification requirements. 
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Table 1 

Definition of the signal regions for each decay mode. 



Mode 




(GeV) 


AM (GcV/c2) 


r~ 


' — > e~7r"'"7r~ 


-0.10 


- +0.04 


-0.014 - 


+0.011 


t" 


■ — > e'^7r~7r~ 


-0.09 


- +0.04 


-0.014 - 


+0.011 


t" 


■ — > ii~tt'^tt~ 


-0.07 


- +0.03 


-0.011 - 


+0.011 


t" 


' — /j^'^tt'tt^ 


-0.07 


- +0.03 


-0.011 - 


+0.011 


t" 


' — e"7r''"if ~ 


-0.10 


- +0.04 


-0.013 - 


+0.011 


T~ 


■ — >■ e~7r~if''" 


-0.10 


- +0.04 


-0.012 - 


+0.010 






-0.10 


- +0.04 


-0.014 - 


+0.010 


T~ 


■ ^ e'K+K- 


-0.10 


- +0.04 


-0.010 - 


+0.008 


T" 


■ e+K-R- 


-0.10 


- +0.04 


-0.013 - 


+0.009 


T" 


— ^ Ijl~''k'^K~ 


-0.08 


- +0.03 


-0.009 - 


+0.009 


t" 


— > ^"tt^K'^ 


-0.07 


- +0.03 


-0.009 - 


+0.009 


T~ 


— ^ IJ^'tt~K~ 


-0.08 


- +0.03 


-0.009 - 


+0.009 




- n'K+K- 


-0.07 


- +0.03 


-0.007 - 


+0.008 


r~ 


H+K-R- 


-0.07 


- +0.03 


-0.007 - 


+0.007 


r~ 




-0.10 


- +0.04 


-0.015 - 


+0.012 


r~ 


■ ^ e-K*(892)0 


-0.10 


- +0.04 


-0.013 - 


+0.011 




■ ^ e"^* (892)0 


-0.08 


- +0.04 


-0.012 - 


+0.010 




" — ^ e"^ 


-0.09 


- +0.03 


-0.010 - 


+0.008 


T~ 




-0.07 


- +0.03 


-0.011 - 


+0.011 


T" 




-0.08 


- +0.03 


-0.009 - 


+0.010 


T~ 


■ ^ /x-^* (892)0 


-0.08 


- +0.03 


-0.009 - 


+0.009 


T~ 




-0.08 


- +0.03 


-0.007 - 


+0.007 



We verify that this shape does not change when additional selection criteria 
are applied by comparing the distribution before all selections and after each 
of them is successively imposed. It is known that K/n separation is different 
for data and MC. We measure identification efficiency and fake rate for both of 
them using cahbration samples obtained from D* — > D^[K ~7r'^]7r decays and 
apply a correction to the MC distributions. The AM distributions of data and 
expected background for the r~ e~Ti^Ti~ and r~ fj~Tr~^7r~ modes after 
the corrections are shown in Figure 3. The number of expected background 
events in the signal region for each studied decay mode is given in the fifth col- 
umn of Table 2. The uncertainty in the background expectation is dominated 
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Fig. 3. AM distributions of experimental data (points with error bars) and expected 
background (histograms) for a) r — > e~7r+7r^ and b) r — > /i^vr+vr^ modes after 
all selection criteria except for AM. Particle identification corrections (see text) 
have been applied to all MC distributions. Different patterns of the background 
histograms correspond to various kinds of the background: uds continuum (dark), 
r-pair generic decays (shaded) and two-photon process (blank). MC histograms are 
cumulative and show general agreement between the data and MC. 

by MC statistics. 

The AE* vs. AM plots for the experimental data for all decay modes are 
shown in Figure 4. The numbers of events observed in the signal regions are 
listed in the sixth column of Table 2. They are consistent with those expected 
from background distributions. We set the upper limits Sq on the number of the 
signal events at 90% CL using the prescription of Feldman and Cousins [18]. 
The main systematic uncertainties on the detection efficiency come from track 
reconstruction (1.0% per track), electron identification (1.1% per electron), 
muon identification (5.4% per muon), kaon/pion separation (1.0% per kaon 
and pion), trigger efficiency (1.4%), statistics of the signal MC (1.0%) and 
uncertainties of the branching fractions of vector meson decays (1.2% for (p 
K'^K~). The uncertainty on the number of r-pair events mainly comes from 
the luminosity measurement (1.4%). The systematic uncertainties due to the 
angular distribution of LEV r decays are summarized in the second column 
of Table 2, together with the total systematic errors Ae/e, which are used to 
evaluate the total errors of the sensitivities. To evaluate Sq, we use the Poisson 
probability density function, which is convolved with the uncertainties on the 
sensitivities and expected background assuming the Gaussian shape [19,20]. 

Upper limits on the branching fractions B are calculated as B{t — ^ ihh) < 

° and B{t^ £V°) < — — — ° , where N^-j. is the total number of 
zNtt^Bi zN-j-T-eBiBy 

the r-pairs produced, Bi is the inclusive branching fraction for the 1-prong 
decay of r and By is the branching fraction of vector meson decay into charged 
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Fig. 4. IS.E* vs. AM experimental data distribution after all selection criteria. The 
boundaries of the signal region are illustrated by solid boxes. 
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Table 2 

Summary of detection efficiency e and its systematic error Ae/e, background ex- 
pectation, number of observed events and 90% C.L. upper limits on the branching 
fractions (BFs). For systematic uncertainty, LFV stands for the error due to the 
angular distribution of LFV decays, while "Total" combines LFV and all other 
errors. 



Mode 


Ae/e 


(%) 


Detection 


Expected 


Observed 


Upper limit on 




LFV 


Total 


efficiency e (%) 


background 


events 


BF (90% CL) 


T~ — > e~TT'^'IT~ 


5.3 


7.5 


5.30 


2.62±1.07 


6 


7.3x10"''' 


T~ — > e"'~7r~7r~ 


2.3 


5.8 


5.14 


0.00±0.26 


1 


2.0x10"''' 


T~ — > /x~7r"^7r~ 


2.1 


8.8 


4.37 


0.76±0.26 


2 


4.8x10"''' 


T~ — > fi'^Tr~Tr~ 


7.7 


11.5 


4.44 


0.73±0.30 


1 


3.4x10"''' 


T~ — > e'~'K'^ K" 


20.5 


21.2 


3.99 


0.91±0.25 


3 


7.2x10"^ 


T~ — > e^TT^ 


17.4 


18.2 


4.11 


1.27±0.41 





1.6x10"^ 


T~ e~^Tr~ K~ 


12.8 


13.9 


4.03 


0.74±0.22 





1.9x10"^ 


T- e-K-K+ 


21.9 


22.5 


3.12 


0.34zb0.20 





3.0x10"^ 


T" e+R-R- 


5.4 


7.6 


3.06 


0.09±0.07 





3.1x10"^ 


T~ — ^ ii'~Tr'^R~ 


15.8 


18.0 


3.43 


2.35±0.44 


1 


2.7x10"^ 


T~ — > llTlT^ R'^ 


19.1 


20.9 


3.32 


1.85±0.32 


3 


7.3x10"^ 


T~ — > fi^Tr~R~ 


25.4 


26.8 


3.53 


2.53±0.38 


1 


2.9x10"^ 


r~ I1-R-R+ 


8.7 


12.2 


2.76 


0.48±0.19 


2 


8.0x10"^ 


T- n+R-R- 


38.2 


39.2 


2.70 


0.09±0.06 





4.4x10"^ 


T~ — >■ e~p^ 


5.3 


7.5 


5.03 


2.55±1.04 


5 


6.5x10"^ 


T- e-R*{m2f 


17.4 


18.2 


4.12 


0.76±0.34 





3.0x10"^ 


r- ^ e-^* (892)0 


20.5 


21.2 


3.68 


0.16±0.10 





4.0x10"^ 


t" — > e~(j) 


21.9 


22.5 


2.94 


0.04±0.04 





7.3x10"^ 


r- /x-pO 


2.1 


8.8 


4.40 


0.26±0.12 





2.0x10"^ 


r- ^ /i-if* (892)0 


19.1 


20.9 


3.61 


0.37±0.14 





3.9x10"^ 


r- ^ /x-X* (892)0 


15.8 


18.0 


3.42 


0.49±0.19 





4.0x10"^ 




8.7 


12.2 


2.68 


0.00±0.18 





7.7x10"^ 



hadrons. We use Bi = (85.35 ± 0.07)% and = (49.2 ± 0.6)% from the 2005 
web update of Ref. [13]. The resulting upper limits on the branching fractions 
are summarized in the last column of Table 2. 
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4 Discussion 



Our final results for 90% upper limits on the branching fractions for the r — > 
ihh modes, shown in Table 2, are in the range (1.6 — 8.0) x 10~^. Our results 

for r — > ihh and r iV^ modes arc one order of magnitude more restrictive 
than those obtained in the CLEO experiment [4] except for a few decay modes. 
The results for r~ — > ihh modes are comparable to those from BaBar [5]. 

While in many extensions of the Standard Model the predicted values of the 
branching fractions of LFV decays are very small and out of reach for current 
experiments, some models allow enhancements of such decays to a level very 
close to the experimentally accessible range at the B factories. For example, 
detailed analysis of various LFV r lepton decays in the framework of MSSM 
show that at small tan/3 and appropriate values of other model parameters 
the branching ratio of the t~ — > ii~p^ decay can be as high as 10~^ [21]. 

The improved sensitivity to rare r lepton decays achieved in this work can be 
used to constrain the parameters of models with heavy Dirac neutrinos [3]. 
In this model the expected branching fractions of various LFV decays are 
evaluated in terms of combinations of the model parameters. These combi- 
nations, denoted y^^ and for r decays involving an electron and a muon, 
respectively, can vary from to 1. Our best 90% upper limit for the modes 
with an electron, y^^ < 0.24, can be set from the r~ — > e~p° decay. The best 
corresponding limit for modes with a muon studied in this work is set from 
the T~ fJ'~P^ mode, < 0.38. This bound is more restrictive than any 
other limits set on LFV decays of the r. 

Our results can also be used to constrain energy scale of new physics in models 
with dimension-six effective fermionic operators that induce r — n mixing [22]. 
From our upper limits for the branching fractions of r~ t~ fi~(f) 

and ^ p-K*{892f the bounds Auu,dd > 29.4 TeV, A,, > 24.8 TeV and 
Ads > 26.8 TeV, respectively, can be obtained for the models with vector 
operators. 



5 Summciry 

We have searched for LFV decays r — > ihh and r — > iV^ using a 158.0 fb~^ 
(140.9 X 10^ T-pair events) data sample in the Belle experiment. No evidence 
for a signal of these decay modes is observed and upper limits on the branch- 
ing fractions are set in the range (1.6 — 8.0) x 10^^, which are one order of 
magnitude more restrictive than those previously obtained by the CLEO ex- 
periment. For T — > ihh modes, the results are comparable to recent limits 
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from BaBar [5]. 
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